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Abstract 

We present highly robust, optimal control-based shaped pulses designed to replace all 90° and 180° hard 
pulses in a given pulse sequence for improved performance. Special attention was devoted to ensuring that the 
pulses can be simply substituted in a one-to-one fashion for the original hard pulses without any additional 
modification of the existing sequence. The set of four pulses for each nucleus therefore consists of 90° and 
180° point-to-point (PP) and universal rotation (UR) pulses of identical duration. These 1 ms pulses provide 
uniform performance over resonance offsets of 20 kHz ( 1 H) and 35 kHz ( 13 C) and tolerate reasonably large 
radio frequency (RF) inhomogeneity/miscalibration of ±15% ( 1 H) and ±10% ( 13 C), making them especially 
suitable for NMR of small-to-medium-sized molecules (for which relaxation effects during the pulse are 
negligible) at an accessible and widely utilized spectrometer field strength of 600 MHz. The experimental 
performance of conventional hard-pulse sequences is shown to be greatly improved by incorporating the new 
pulses, each set referred to as the Fantastic Four (Fanta4). 
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1. Introduction 

Hard rectangular pulses are the fundamental 
element of all multi-dimensional NMR pulse se- 
quences. They are simple to use and extremely 
versatile. The transformations required for multi- 
dimensional spectroscopy — excitation, flip-back, in- 
version, and refocusing — are all possible using this 
simplest of pulses. Moreover, the performance of 
a single hard pulse is approximately ideal over a 
reasonable range of resonance offsets (proportional 
to the RF amplitude of the pulse) and variation in 
RF homogeneity (±10%) relevant for high resolu- 
tion spectroscopy. 

However, there is considerable room for improv- 
ing pulse sequence performance. RF amplitude 
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is limited in practice and cannot be increased to 
match the increased pulse bandwidth needed at 
higher field strengths. The problem of large chem- 
ical shift can be solved by dividing the spectral 
region and performing multiple experiments. But 
this is time consuming, and unstable samples can 
create additional problems. Even at lower field 
strengths, relatively small errors produced by a sin- 
gle pulse can accumulate significantly in multipulsc 
sequences. 

The use of shaped pulses that address particular 
limitations of hard pulses can improve performance 
[D-U8]. but complex multipulse sequences are mas- 
terpieces of timing and synchronized spin-state evo- 
lution. Replacing e.g. a given pair of rectangular 
pulses (c.f. Fig. [iji) in a sequence with a better 
performing pair of shaped pulses with different du- 
rations (c.f. Fig. [TJd) typically requires a nontrivial 
redesign of the pulse sequence. Additional pulses 
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and delays are required to reestablish the timing 
and refocusing that achieve the goals of the orig- 
inal hard-pulse sequence. In order to avoid these 
complications, here we developed a set of shaped 
pulses with identical durations, such that the re- 
placement of rectangular pulses by shaped pulses 
does not introduce additional delays (c.f. Fig. [lp). 

The goal of the present work is to provide a 
fundamental set of better-performing pulses that 
can simply replace, in one-to-one fashion, all the 
hard 90° and 180° pulses in any existing NMR se- 
quence. This includes important sequences such as 
HSQC, HMBC, HMQC, INADEQUATE, COSY, 
and NOESY, which are some of the most basic pulse 
sequences used for finding correlations between nu- 
clei within a given molecule. Since hard pulses 
perform a universal rotation (UR) about a given 
fixed axis for any orientation of the initial magne- 
tization, better-performing UR 90° and UR 180° 
pulses are, in principle, all that are needed. How- 
ever, UR pulses are only strictly necessary for refo- 
cusing or combined excitation/flip-back. Point-to- 
point (PP) pulses, which transform only one specific 
initial state to a desired final state, are sufficient for 
excitation, flip-back, or inversion, and these can be 
designed more easily and with better performance 
for their specific task than UR pulses [T9Tf2"l"] . 

We used the optimal control-based GRAPE al- 
gorithm [22] to design 90° and 180° pulses for 
UR and PP transformations applicable to 1 H and 
13 C spectroscopy at 600 MHz, which is currently a 
generally accessible and widely utilized spectrom- 
eter field strength. Pulse bandwidths are there- 
fore 20 kHz for X H (33.33 ppm) and 35 kHz for 
13 C (233.33 ppm). In order to universally imple- 
ment these shaped pulses on any available probe- 
head, the maximum RF amplitude for 13 C spins 
was set to RF m &x = 10 kHz, robust to ±10% 
RF inhomogeneity/miscalibration. For *H nuclei, 
flFmax = 18 kHz was allowed, with tolerance to 
±15% variations. 

A pulse length T p = 1 ms is the minimum length 
necessary to achieve suitably uniform performance 
for the 180° UR transformation over resonance off- 
sets of 35 kHz with i?F m ax limited to 10 kHz. Al- 
though the same performance can be obtained using 
shorter pulse lengths for the other pulse criteria, all 
pulses were designed using the same 1 ms length. 
Each hard pulse in a sequence can then be replaced 
by the corresponding UR or PP pulse without any 
further modification of the sequence. This pulse 
length is most suitable for spectroscopy of small- 



to medium-sized molecules, where relaxation effects 
are neglibile during the pulse. The final product is 
two sets of four pulses, each set whimsically referred 
to as the Fantastic Four (Fanta4) due to the signif- 
icant improvement they provide in pulse sequence 
performance. 

2. Optimization 

The GRAPE algorithm for pulse optimization is 
discussed in detail in 122 . Further details of its 
application can be found in the cited references on 
optimal control [THl |2"3TI3"T] . and in particular the 
references on UR pulses [T9"H2"1"] . A quality fac- 
tor, $, for pulse performance is defined which, in 
turn, enables an efficiently calculated gradient for 
iterative improvement of pulse performance. Most 
generally, the quality factor is a quantitative com- 
parison between the state of the system and some 
desired target state. The gradient therefore also de- 
pends on the system state which, for a single spin, 
is given by the magnetization M. Modifications in 
the basic algorithm that are required for two-spin 
systems requires some elaboration. 

2.1. Two-spin systems 

Evolution due to heteronuclear J-coupling during 
a pulse must be considered during the optimization 
of the Fanta4 pulses. At the simplest level, a refo- 
cused pulse with no chemical shift evolution during 
the pulse produces no heteronuclear J-coupling evo- 
lution when applied to a single spin [29] . However, 
this is not the case when pulses are applied simul- 
taneously to X H and 13 C, where Hartmann-Hahn 
transfer must be considered. In principle, this is 
readily addressed by optimizing shaped pulses for a 
coupled two spin-1/2 system simultaneously. 

However, optimization of shaped pulses for cou- 
pled two-spin systems over a range of resonance off- 
sets and tolerance to RF inhomogeneity is compu- 
tationally very expensive. For example, a one-spin 
system with an offset range of 20-35 kHz digitized 
in 250-500 Hz increments requires on the order of 
100 offsets (choosing rounded numbers for simple il- 
lustration). Tolerance to RF inhomgeneity of ±10- 
15% can be included using ~ 10 RF increments, 
giving an order-of-magnitudc total of 1000 combi- 
nations. The time evolution and cost 4> must be 
calculated for each combination given a particular 
trial pulse in the iterative optimization procedure. 
Adding a second spin, the number of combinations 
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becomes 10 6 , increasing the computation time by 3 
orders of magnitude. The actual time can easily be 
4 orders of magnitude larger than the single-spin 
case when all actual additional factors of 2-3 are 
included. Even parallel processing (as we have uti- 
lized for a decade) requires a nontrivial amount of 
processing time. Details of the coupled two-spin- 
1/2 optimizations and the resulting pulses are pro- 
vided in the Appendix (see section 7.2 1. 



The optimization algorithm begins with a trial 
pulse. In the absence of a known pulse that per- 
forms reasonably well as a candidate for the op- 
timization, the initial pulse is typically generated 
as a sequence of random amplitudes and phases. 
Different initial pulses most generally result in dif- 
ferent final pulses which can vary in performance, 
so some experimentation is necessary using multi- 
ple trial pulses to determine the ultimate quality 
factor that is physically attainable for a given set 
of performance criteria. 

2.2. Single- spin proxy 

The physical performance limits of pulses applied 
to single-spin systems have been studied in detail 
for broadband excitation and inversion, both for 
fixed RF max [55] and fixed RF power [3U]. A com- 
plete and thorough characterization of UR pulse 
performance limits is also available [32[ [33] . 

We used these single-spin results as a guide for 
what might be possible for two-spin systems. As an 
alternative to multiple, lengthy two-spin optimiza- 
tions, we adopted an intuitive approach based on 
the likelihood that pulses separately and indepen- 
dently optimized for 1 H and 13 C (with different rf 
amplitudes) using a single-spin protocol would not 
satisfy the frequency (or phase modulation) match- 
ing conditions for sufficient Hartmann-Hahn trans- 
fer. 

We therefore also approximated the coupled 
two-spin-1/2 problem as two independent, non- 
interacting single spin-1/2 problems, each of which 
demands less computation by several orders of mag- 
nitude compared to the two-spin case. We then 
simulated simultaneous application of the resulting 
pulses to 1 H and 13 C using the full two-spin-1/2 
system with a J coupling of 200 Hz. 

For example, a PP 90° shaped pulse along the 
x-axis should rotate initial z-magnetization to the 
— y-axis. In the case of a coupled two-spin system, 
a combination of a PP 90° pulse on 1 H (/-spin) 
and 13 C (iS-spin) should rotate initial I z and S z to 
I y and S y with maximum fidelity and a minimum 



of unwanted terms such as 2I y S z or 2I z S y . The 
best shaped pulses which fell within the threshold 
value of fidelity were included in the Fanta4 pulse 
set. Further details of Fanta4 pulse selection using 
this protocol are provided in the Appendix (section 
7Tf. 

We also note that if only one pulse is applied 
to either of the spins, transverse magnetization of 
the other spin can evolve under the chemical shift 
Hamiltonian. We therefore require a "Do-Nothing" 
pulse which rotates magnetization by an arbitrary 
multiple of 360° , which can be realized by two iden- 
tical PP inversion pulses with a duration of 500 /is 
each. 



3. Experiments 

We chose the HSQC (Fig. [2) and HMBC (Fig. [3) 
experiments to demonstrate the implementation 
and performance advantages of Fanta4 pulses. 
HSQC is widely used for recording one-bond corre- 
lation spectra between two heteronuclei. HMBC is 
mostly used for correlating heteronuclei connected 
by multiple bonds, mostly 2-4 bonds. Hard pulses 
were replaced by corresponding Fanta4 pulses in 
these pulse sequences. For the most efficient re- 
placement, the task of each hard pulse in Figs. [2^, 
and [3^,, is analyzed and classified according to the 
following four pulse types: UR 90°, PP 90°, UR 
180° and PP 180° (represented in Figs. |]b and ^> 
by solid black bars, black bars with white horizon- 
tal stripes, white bars and white bars with black 
horizontal stripes, respectively). For example, the 
first *H 90°2;-pulse acts on initial z magnetization 
and hence can be replaced by a PP 90° pulse that 
transfers I z to —I y . In the HSQC experiment, 
the following 1 H-180°pulse is a refocussing element 
that needs to be implemented as a UR 180°pulse, 
whereas the simultaneously applied 13 C-180°pulse 
only needs to invert 13 C spins and can be replaced 
by a PP 180°pulse. Note that in the back transfer 
with sensitivity enhancement, some of the 90° pulses 
(represented by solid bars in Fig. [2]b) also need to 
be implemented as UR 90°pulses as they simultane- 
ously need to act on two orthogonal spin operators. 

All experiments were performed at 298° K 
on a Bruker 600 MHz AVANCE III spectromter 
equipped with SGU units for RF control and lin- 
earized amplifiers, utilizing a triple-resonance TXI 
probehead and gradients along the z-axis. Further 
details of individual pulse shapes and performance 
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as a function of resonance offset are provided in the 



Appendix (see section 7.1.1 ) 



4. Results and Discussion 

4-1- HSQC: Testing with Sodium Formate 

13 C-labeled Sodium formate (Fig. |4^), forming 
a simple two-spin system with one proton and 
one carbon spin (J = 197 Hz) dissolved in D 2 
was used to compare the performance of conven- 
tional, and adiabatic-HSQC [53H3"6] (where 13 C 
180° hard pulses were replaced by adiabatic inver- 
sion and composite refocusing pulses constructed 
using Chirp60 [TH [15] from Bruker library [37] ) 
to the Fanta4-HSQC sequence at different chemi- 
cal shifts and RF miscalibrations. To reduce the 
effects of RF field inhomogencity, a small volume of 
approximately 40 fA of sample solution was placed 
in a 5 mm Shigemi limited volume tube. A proton- 
excited and detected HSQC (coupled during acqui- 
sition) gives a proton doublet with respect to car- 
bon. 

1D-HSQC experiments were performed incre- 
menting the carrier frequency successively over the 
offset range in the 13 C dimension. The correctly 
calibrated RF amplitude is applied to X H and 13 C 
in the comparison between conventional, adiabatic, 
and Fanta4 pulses shown in Figs. [5^, [5p, and [Hp. 
The conventional HSQC produces signal intensity 
within 80% of the on-resonance signal of the cali- 
brated pulse over at most an offset range of ±5 kHz, 
far short of the desired 35 kHz bandwidth. RF 
pulse errors accumulate during the sequence, result- 
ing in degraded signal intensity with considerable 
phase errors, in contrast to the nearly uniform per- 
formance of HSQC sequence with adiabatic 180° 
pulses on 13 C and the Fanta4 implentation over 
the 35 kHz offset range. The small loss of signal 
(~ 10%) on resonance for the Fanta4 HSQC com- 
pared to the hard-pulse HSQC is more than com- 
pensated by signal gains and more uniform perfor- 
mance over the desired bandwidth using Fanta4. 
The contrast is especially striking given that the 
peak amplitudes of the hard pulses were higher than 
the Fanta4 peak amplitudes by almost a factor of 
two for 13 C and a factor of 1.6 for : H. 

For the case of adiabatic 180° pulses, peak ampli- 
tude was 10 kHz with durations of 2 ms for refocus- 
ing and 0.5 ms for inversion pulses. The adiabatic 
HSQC was included to show the performance that is 
available by adjusting pulses for a specific applica- 
tion rather than using the Fanta4 pulses, which are 



designed to be applied more generically. As noted 
in the Introduction, PP pulses can be designed for 
a specific task and performance criterion that are 
shorter than the corresponding UR pulses, but re- 
quire additional adjustment to the pulse sequence. 
In the present comparison, the adiabatic refocusing 
pulse applied to 13 C is twice as long as the Fanta4 
UR 180° pulse, resulting in the better performance 
shown in Fig. pp. However, miscalibrating the RF 
amplitude on 1 U and 13 C by -15% and -10%, 
respectively, in the conventional/adiabatic/Fanta4 
comparison of Figs. [5p>, [5pl and[5p only marginally 
affects the Fanta4 implementation, while the con- 
ventional HSQC is almost unusable, and adiabatic 
HSQC performed with ~ 50% loss in signal inten- 
sity One could further optimize the performance of 
the HSQC sequence by optimizing all the 1 H and 
13 C pulses and pulse timings, which is not the focus 
of the current article. 

4-2. HSQC: Testing with Hydroxycitronellal 

We also implemented the 

conventional/adiabatic/Fanta4-HSQC sequences 
on the more complex molecule Hydroxycitronellal 
(Fig. lib). (In contrast to the experiments with 
13 C-labeled Sodium formate with a Shigemi tube, 
a normal 5 mm NMR tube was used with a sample 
volume of about 600 /A.) The molecule consists 
of a long chain of hydrocarbons with a hydroxyl 
group on one end and an aldehyde moiety at other 
end. Hydrocarbons resonate around 19 ppm, and 
a carbon bonded to oxygen (aldehyde moiety) 
resonates at 202 ppm, which corresponds to a 
total 13 C offset range of 27.3 kHz on a 600 MHz 
spectrometer. 

Even with correctly calibrated RF amplitude, the 
conventional HSQC experiment is unable to excite 
the 13 C nuclei at very large offsets and results in 
poor signal-to-noise ratio in Fig^Bp, compared to the 
adiabatic and the Fanta4 HSQC shown in Fig. [6b 
and Fig. [6p respectively. The limitations of hard- 
pulse and adiabatic 180° pulse implementations of 
HSQC are further emphasized in the comparison 
between Fig.[6ji, Fig.[6p and Fig. [Bf, with the RF 
amplitude miscalibrated on both nuclei as in the 
previous section. 

4-3. HMBC: Testing on a complex molecule 

A Furan with substrate l-methoxy-4- 
methylbenzene (4-(5-(4-mcthoxyphcnyl)-3- 
methylfuran-2-yl)butanal) (Fig. |4p), synthesized 
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by a domino reaction (see table 3, entry 2 in [3"5]). 
was used to test the HMBC pulse sequence. 

It is a relatively small molecule with a large num- 
ber of long-range J couplings between 1 H and 13 C 
spins. The required 13 C bandwidth is 28.8 kHz at 
600 MHz. The HMBC pulse sequence with hard 
pulses (red) shows poor performance at large offsets 
compared to the Fanta4-HMBC (black) (Fig. [7|). 
Traces from this 2D HMBC spectrum are compared 
in Fig. [HJ providing further detail of the signal-to- 
noise enhancement available using Fanta4, showing 
gains of up to a factor of three in signal-to-noise 
ratio. 



5. Conclusion 



7. Appendix 

As discussed in the main text, the sets of Fanta4 
pulses presented for *H and 13 C spins were opti- 
mized assuming single, uncoupled spins. In sec- 
tion [7tJ we first explain the post-selection protocol 
that was used to conduct a thorough search and 
select the best combination of pulses (in terms of 
their performance for coupled spins) that were in- 
dividually optimized for uncoupled spins. In sec- 
tion |7.2[ we discuss the far more time consuming 
approach of direct pulse optimization for coupled 
spins where spin-spin couplings are explicitly taken 
into account. Due to its extreme computation de- 
mands, this procedure was only used to develop a 
representative pulse. 



We have designed two sets of shaped pulses, re- 
ferred to as Fanta4 pulses, that can be used to re- 
place all 90° and 180° hard pulses on *H and 13 C, 
respectively, in conventional pulse sequences for im- 
proved performance. HSQC and HMBC experi- 
ments were provided to show that each hard pulse 
in a sequence can be easily replaced with the corre- 
sponding Fanta4 shaped pulse without further mod- 
ifying the existing pulse sequence. Compared to 
rectangular pulses, Fanta4 pulses provide far more 
robust performance with respect to large frequency 
offsets, RF inhomogeneity and/or RF miscalibra- 
tion. The duration of the current generation of 
Fanta4 pulses is 1 ms, which renders them suitable 
for small and medium sized molecules with moder- 
ate relaxation values. 



6. Outlook 

Shorter pulses: Current Fanta4 pulses are 1ms 
long, making them most suitable for spectroscopy 
of small-to medium-sized molecules, where relax- 
ation effects are neglibile during the pulse. For ap- 
plications to large molecules with considerable re- 
laxation effects, we are in the process of designing 
shorter pulses. 

Two-spin optimizations: The two-spin optimiza- 



tion process detailed in the Appendix (7.2) showed 
promising results. To address the previously 
noted problem of significantly increased computa- 
tion time, we will consider a combination of a new 
optimization algorithm(s) and parallel programing. 



7.1. Fanta4 pulse selection 

Here we explain the steps of the post-selection 
protocol that we used to select from a large set of 1 H 
and 13 C PP and UR pulses that were found based 
on single-spin optimizations to obtain the final set 
of Fanta4 pulses presented in this paper. Finding 
the best combination of pulses based on their sim- 
ulated performance in coupled two-spin systems is 
a nontrivial combinatorial problem which becomes 
tractable using the following procedure. 

1. From a pool of hundreds of optimized 1 H 
and 13 C PP 90°, PP 180°, UR 90°, and UR 180° 
pulses (optimized starting from different random se- 
quences) with a duration of 1 ms, between four and 
13 individual pulses with fidelity above 0.9999 for 
uncoupled spins were chosen as final candidates for 
each class of *H and 13 C PP 90°, PP 180°, UR 90°, 
and UR 180° pulses. The set of final candidate 
pulses is shown schematically in Fig. [9^, (assuming 
for simplicity four pulses in each class). 

2. For each pair of candidate X H and 13 C pulses, 
the evolution of the fifteen orthogonal initial den- 
sity operators p^'(0) representing a two-spin-1/2 
system (given by the Cartesian basis operators 
p«(0) = I x , p( 2 >(0) = I y , p^(0) = I„ p^(Q) = 
S x , p^(0) = S y , p^(0) = S z , pW(0) = 2I X S X , 
p (s \0) = 2I y S x , ••■ p< 15 )(0) = 2I Z S Z ) was simu- 
lated assuming a J-coupling of 200 Hz and a nom- 
inal RF amplitude of 18 kHz and 10 kHz for *H 
and 13 C pulses, respectively. This is repeated for 
all combinations of 81 offsets v H and 141 offsets v 
resulting from a digitization of the 1 H offset range 
of 20 kHz and the 13 C offset range of 35 kHz in 
steps of 250 Hz. 
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3. For each combination of the four 1 H and 13 C 
classes, we partitioned the resulting 15 Cartesian 
product operators for a two-spin system in three 
groups: desired terms (D), undesired terms (U) 
and terms which can be ignored (I). As the indi- 
vidual pulses have excellent performance for single 
spins, special attention is given here to potential ad- 
verse effects due to the J coupling, such as (partial) 
Hartmann-Hahn transfer in the two-spin system. 

For example, consider the case of the simulta- 
neous application of a PP 90° 1 H (spin /) pulse 
(transferring I z to —I v for an uncoupled spin /) 
and a PP 90° 13 C (spin S) pulse (transferring S z 
to — S y for an uncoupled spin S) in the presence 
of J coupling. If the initial density operator term 
is e.g. p^(0) — I z , obviously the desired final op- 
erator is —Iy Hence, the expectation value of the 
operator I y is considered to be a desired term (D), 
taking into account the correct phase according to 
the specification of the pulse. 

The expectation values of all other terms should 
be as small as possible and hence are undesired (U). 
However, note that a PP 90° 1 H (spin I) pulse that 
transfers I z to —I y does not necessarily leave I x 
invariant - in fact the improved performance of PP 
90° pulses results from the additional degree of free- 
dom related to the fact that the final term can be 
anywhere in the x-z plane. Hence, if the initial den- 
sity operator term is pW(0) = I x , the expectation 
values of both I x and I z are arbitrary and will not 
play a role in the pulse sequence, otherwise a UR 
pulse should have been chosen at this point in the 
sequence. Therefore, for p^'(0) = I x , the expecta- 
tion values of the operators I x and I z can be ignored 
(I) , whereas all remaining operators are undesirable 
(U) as the initial I x operator should not leak to 
any of the remaining 13 product operator terms as 
a result of coupling evolution. 

For all combinations of X H and 13 C pulse classes, 
the tables shown in Figs . |l0| - [13] summarize the de- 
sired terms (D), the undesired terms (U), and the 
terms of the final density operator that can be ig- 
nored (I) for all considered initial Cartesian product 
operators p (1) (0), ■•■ p (15) (0). 

With this classification of expectation values into 
terms that are desired, undesired or that can be ig- 
nored, the following quality factor is calculated for 
each pair of the candidate 2 H and 13 C pulses to 
quantify their combined performance in the pres- 
ence of the J coupling: 



where Q\j and are the individual quality fac- 
tors for the desired (D) and undesired (U) terms for 
this pulse pair 

and and are the averages of and 

over all combinations of offsets v H and v c . 

The quality factor Q\f is defined as as the aver- 
age of the absolute values of the expectation values 
for the desired transformations given in the corre- 
sponding table. For example, if the candidate pulse 
i is a PP 90° for *H and pulse j is a PP 90° for 13 C, 



the corresponding table is found in Fig. 10 (top left) 
and in this specific case 



Q 
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(\(Iy) {3) \ 
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where (A)<p> = Tr{A*pW(T)} with p^{0) = I z 
as defined above in step 2. Similarly, the quality 

factor Q^p is defined as the average of the absolute 
values of the expectation values for the undesired 
transformations given in the corresponding table. 

Using Eq. ([!]), we can calculate the quality fac- 
tors Qij for all pulse pairs ij and fill the master 
table represented in Fig. [9^. 

4. In this step, all possible sets of Fanta4 pulses 
(consisting of one PP 90°, one PP 180°, one UR 
90°, and one UR 180°for both 1 H and 13 C) are con- 
structed and the overall performance of each set is 
quantified (cf.Fig.J9j3). The quality factor for each 
set is simply given by the average of the quality 
factors Qij of all pulse pairs in the set, where the 
values of Qij have been calculated in step 3 and can 
be taken from the master table (Fig. [9^,). Finally, 
the set with the best overall quality factor was cho- 
sen as the set of Fanta4 pulses presented in this 
paper. The individual quality factors Qij for the 
best set of Fanta4 pulses are summarized in Table 
1. 

The set of best pulses reduces the effect of het- 
eronuclear J coupling during the shaped pulse. 

7.1.1. Fanta4 pulse shapes and excitation profiles 

All experiments were implemented on a Bruker 
600 MHz AVANCE III spectromter equipped with 
SGU units for RF control and linearized amplifiers, 
utilizing a triple-resonance TXI probehead and gra- 
dients along the z-axis. Measurements are the resid- 
ual HDO signal using a sample of 99.96% D2O 
doped with CUSO4 to a T\ relaxation time of 100 ms 



Qij Q 



(D) 

ij 



Q 



(U) 



(1) 



at 298° K. For the X H pulses shown in Figure 14 



signals are obtained for offsets between -11.1 kHz 
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Table 1: Q from Eq. [ T| is summarized for the current set of 1 H and 13 C Fanta4 p ulses. 



PP H 90° 
PP H 180° 

UR H 90° 
UR H 180° 



PP C 90° PP C 180° UR C 90° UR C 180° 



0.9603 
0.9796 
0.9832 
0.9786 



0.9762 
0.9767 
0.9753 
0.9753 



0.9725 
0.9717 
0.9710 
0.9657 



0.9679 
0.9701 
0.9634 
0.9665 



to 11.1 kHz in steps of 200 Hz at an ideal RF am- 
plitude with RF max of 18 kHz (Fig. p). For the 



13 C pulses shown in Figure 16 signals are obtained 
at offsets between -18.5 kHz to 18.5 kHz in steps 

of 



of 200 Hz at an ideal RF amplitude with RF n 
10 kHz (Fig. 



171. To reduce the effects of RF field 



inhomogeneity, approximately 40 p\ of sample solu- 
tion was placed in a 5 mm Shigemi limited volume 
tube. The duration of each Fanta4 pulse is 1 ms. 

7.2. Optimization of pulse pairs for two coupled 
heteronuclear spins 1/2 

In the previous section, the protocol was ex- 
plained that was necessary to select the set of 
Fanta4 pulses presented in this paper based on indi- 
vidual pulses optimized for uncoupled spins. Here, 
we show how the optimal control approach can be 
used to directly optimize pulse pairs for two J- 
coupled hetero-nuclear spin-1/2 nuclei. However, 
as this algorithm is computationally demanding, we 
demonstrate this approach only for a single pulse 
pair. 

7.2.1. Transfer from single initial to single final 
state 

The goal is to find a RF pulse which steers the 
trajectory from a given initial product operator 
term p(0) — 2I a Sp at time t — to a desired tar- 
get state F = 2I 7 Ss at t = T, where a, [3, 7, and S 
correspond to the x, y, and z components of mag- 
netization. This is accomplished by optimizing a 
suitably chosen cost function or performance index, 
$, as discussed further below. 

The state of the spin system at time t is char- 
acterized by the density operator p(t). The equa- 
tion of motion in the absence of relaxation is the 
Liouville-von Neuman equation, 



p(t) = -i[(fl + X]«fc(t)fffc),P(*)] > 



(3) 



k=l 



where Hq is the free evolution Hamiltonian, 
are the RF Hamiltonians corresponding to available 



control fields, and u(t) = (ui(i), Ua(i), u m {t)) 
represents the vector of RF amplitudes, referred to 
as the control vector. For the case of two spins con- 
sidered here, m — 4 for the x- and y-components 
(or, alternatively, amplitude and phase) of the RF 
fields applied to the two different spins. The prob- 
lem is to find the optimal control amplitudes Ukit) 
that steer a given initial density operator p(0) in a 
specified time T to a desired target operator F with 
maximum overlap. 

For Hermitian operators p(0) and F, this overlap 
can be measured by the standard inner product 



(F\p(T)) = Tr{F\p(T)}, 



(4) 



where the operator Tr returns the trace (sum of di- 
agonal elements) of its argument. Hence, the per- 
formance index, $, of the transfer process can be 
defined as 



$ = (F\p(T)). 



(5) 



For the full treatment of the optimization procedure 
we refer to Ref. |39j . 

Digitizing the pulse in N equal steps indexed by 
j = 1,2,3, ... ,N, the basic GRAPE algorithm for 
this cost is 

1. Guess initial controls Uk(j)- 

2. Starting from po, calculate p{j) = 
Vy U lPo U\ ■■■U] for all j < N. 

3. Starting from \m = F, calculate X(j) = 
U} +1 ■ ■ ■ U^ N FU N ■ ■ ■ U j+1 for all j < N. 

4. Evaluate the gradient 

5$/5u k {3) = -{\ j \iAt[H k ,p j \) 

and update the m x N control amplitudes 
u k (j)- 

5. With these as the new controls, go to step 2. 
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The algorithm is terminated if the change in 
the performance index $ is smaller than a chosen 
threshold value. 

Treatment of non-Hermitian operators and 
relaxation-optimized coherence transfer can be 
found in Refs. [331 HO] • 

1.2.2. Transfer from two initial to two final states 
Consider a RF pulse which simultaneously trans- 
forms two orthogonal initial operators pi(Q) — 
2I a Sp and P2(0) = 2I^Ss to the respective fi- 
nal states F\ — 2I e S^ and F 2 = 2I n Sg, where 
a, j3, 7, 5, e, £, n, and correspond as above to the 
x, y, and z components. Each state of the spin sys- 
tem is characterized by the density operator p\(t) 
and p 2 {t) at time point t. The Liouvillc-von Ncu- 
man equation for each state is given by 



Pn (t) = -i[(H + J2u k (t)H k ),p n (t)}, 



(G) 



fc=i 



where n = 1, 2, • • • , P, labels the states. 

The overall performance index $ can be defined 

as 



(F 1 | Pl (T)) + (F 2 | P2 (T)) 



(7) 



The modified GRAPE algorithm for two states is 

1. Guess initial controls u k (j) 

2. Starting from p nu , calculate p n j = 
Uj---Ux p n0 U\---U] for all j < N and 
n = l,2. 

3. Starting from A n jy = F n , calculate X n j = 
Ut F n U N --- U j+1 for all j < N and 



U 



71= 1,2. 

4. Evaluate the gradient 

5$/6u k (j) = 

(\u\iAtlHk, Plj ]) + (X 2l \iAt[H k ,p 2j }} 
2 

and update the m x N control amplitudes 
u k {j). 

5. With these as the new controls, go to step 2. 



Figure 18 shows the shape of the pulses optimized 
simultaneously using above algorithm for ^ and 
13 C with J coupling of 197 Hz for the following 
transfers. 



MO) 

p 2 (0) 



2I Z S X 
2I Z S V 



Fl 
F2 



~2IyS Z 



(8) 
(9) 



In this example, the maximum pulse amplitudes for 
both 1 H and 13 C pulses are 10 kHz and the 1 H 
and 13 C offset ranges were ±5 kHz and ±20 kHz, 
respectively. Figure [19] shows the simulations for 
the corresponding transfer efficiencies. 
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Figure 1: This schematic shows a pair of hard pulses with flip angles a and fi, and delays t\ and T2 (Fig. (a)) in an any given 
pulse sequence. If only 13 C rectangular pulses are replaced by shaped pulses (or if shaped pulses of different durations are used 
for 13 C and 1 H), the effect of chemical shift and coupling evolution during the extra delays <5i and 82 have to be taken into 
account (Fig. (b)). However, if rectangular pulses are replaced by Fanta4 pulses of identical duration (Fig. (c)), no additional 
evolution periods have to be taken into account. 
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Figure 2: Proton excited and detected HSQC 134113511411 experiment using conventional 90° and 180° hard pulses (conventional- 
HSQC, Fig. (a)) and Fanta4 pulses (Fanta4-HSQC, Fig. (b)). Phases are x with the exception of <J>i = x -x, <J>2 = x x -x -x, <3>3 
= x x -x -x, $4 = y y -y -y, <& ree = x -x -x x. Delays are A = 1/(4J) for CH groups and Ai = 1/(8J) for all multiplicities. S 
and Si are delays for gradients (Gi and G2) including recovery time, 'tr' is a time reverse shaped pulse. Both pulse sequences 
are practically of identical lengths but the offset-compensated and RF robust Fanta4-HSQC provides higher sensitivity. 
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Figure 3: Proton excited and detected HMBC experiment [42] using conventional 90° and 180° hard pulse (conventional-HMBC 
(Fig. (a))) and Fanat4 pulses (Fanta4-HMBC (Fig. (b))). Phases are x with the exception of $i = x -x, <J>2 = x x -x -x, <3>3 = x 
x x x -x -x -x -x, "Iv ec = x -x x -x -x x -x x. Delays are A = 1/2 J for evolution of long range couplings. 8 is delay for gradients 
(Gi and G2) including recovery time. Both pulse sequences are practically of identical lengths but the offset-compensated and 
RF robust Fanta4-HMBC provides higher sensitivity. 
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Figure 4: (a) Sodium formate, (b) Hydroxycitronellal, and (c) 4-(5-(4-methoxyphenyl)-3-methylfuran-2-yl)butanal from |38l 
are used for comparison of conventional HSQC and HMBC pulse sequences, and Fanta4 -HSQC and -HMBC pulses sequences. 
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Figure 5: Chemical shift (Av) excitation profile comparison of the conventional-HSQC (Fig. (a) and (b)), the adiabatic-HSQC 
(Fig. (c) and (d)), and the Fanta4-HSQC (Fig. (e) and (f)) acquired on a 13 C labeled Sodium formate dissolved in D2O. 
HSQC (coupled during acquisition) gives doublet of a proton with respect to a carbon. Figures (a), (c), and (e) compare 
signal intensity at different 13 C resonance offsets with ideal RF amplitude, and Figures (b), (d), and (f) show the effect of 
RF inhomogeneity/miscalibration for 1 H (-15%) and 13 C (-10%) pulses for different 13 C offsets. In conventional-HSQC and 
adiabatic-HSQC, the ideal RF amplitude on X H hard pulses was 29.51 kHz and 13 C hard pulses was 19.53 kHz, while adiabatic 
shaped (inversion and refocusing) pulses on 13 C were with 10 kHz maximum amplitude. In case of Fanta4-HSQC, shaped 
pulses on X H had a maximum amplitude of 18 kHz, while on 13 C 10 kHz. 
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Figure 6: 2D HSQC projection acquired on 13 C natural abundant Hydroxycitronellal, dissolved in deuterated tetrahydrofuran, 
using conventional (Fig. (a) and (d)), adiabatic (Fig. (b) and (e)), and Fanta4 (Fig. (c) and (f)) -HSQC pulse sequences 
(Fig|2JI. Figure (a) shows the loss in signal intensity at edge of the 13 C dimension compared to Fig. (b), and Fig. (c) with 
nominal RF amplitude. In contrast to conventional and adiabatic HSQC, only the Fanta4 HSQC is relatively unaffected by 
RF miscalibration of -15% on 1 H and -10% on 13 C pulses (compare Fig. (d), Fig. (e), and (f)). Amplitude on 1 H hard pulses 
was 32.26 kHz and 13 C hard pulses was 19.53 kHz in conventional and adiabatic-HSQC. The adiabatic shaped (inversion and 
refocusing) pulses on 13 C were with 10 kHz maximum amplitude. While Fanta4 pulses were with 18 kHz and 10 kHz on X H 
and 13 C respectively. Peaks in 1 H and 13 C dimension are shifted for clear distinction. 1024 X 256 data points were acquired 
with corresponding spectral widths of 230.0 ppm ( 13 C) and 10.0 ppm ( J H). The carrier frequency on X H was 5.3 ppm and 
on 13 C was 115.0 ppm. Two transients per increment gave an overall experiment time of ~ 13.5 min for each of the three 
experiments. 
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Figure 7: A magnitude mode 2D HMBC of a 13 C natural abundant 4-(5-(4-methoxyphenyl)-3-methylfuran-2-yl)butanal 
molecule |38) . dissolved in CDCI3, was acquired using conventional-HMBC (red) and Fanta4-HMBC (black) pulse sequences 
(Fig. [3). The hard pulse version shows loss in signal intensity compared to Fanta4-HMBC (contours are slightly shifted in 13 C 
dimension) near <5 13 C = 10 to 30 ppm and at S 13 C = 202 ppm (encircled). 1024 X 256 data points were acquired with 
corresponding spectral widths of 230.0 ppm ( 13 C) and 10.0 ppm ^H). The carrier frequency on 1 H was 5.4 ppm and on 13 C 
dimension was 115.0 ppm. Forty transients per increment gave an overall experiment time of 6 hrs 25 min for each of the two 
experiments. Amplitude on X H hard pulses was 30.79 kHz and 13 C hard pulses was 19.53 kHz in conventional-HMBC, while 
Fanta4 amplitude were 18 kHz and 10 kHz on X H and 13 C, respectively. 
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Figure 8: (Color online) Traces through cross signals at 21.2, 106.6 and 202.2 ppm are compared for conventional-HMBC (red, 
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Figure 10: (Color online) This figure shows the combinations of PPh 90° pulse with rest of 13 C pulses. For simulation of pair 
of pulses on coupled two spins-1/2 system, each component of cartesian product operator is considered as initial state and all 
components are detected. The desired(D, in red), undesired (empty box), and terms which can be ignored (I) are indicated in 
each case and described more fully in section [7. 1| 
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Figure 11: The combinations of PPjj 180° pulse with the 13 C pulses. Compare Figure 
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ure 12: The combinations of URh 90° pulse with the 13 C pulses. Compare Figure [To] 
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gure 13: The combinations of URjj 180° pulse with the 13 C pulses. Compare Figure 
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Figure 14: The amplitude and phase of 1 H pulses are plotted as a function of the pulse duration. All are 1 ms long with a 
maximum amplitude of 18 kHz. 
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Figure 15: The experimental performance of 1 H pulses is plotted as function of the offset Av at ideal RF amplitude with a 
RFmnx of 18 kHz. Experiments were aquired on residual proton in HDO for an offset range of 22.2 kHz in steps of 200 Hz. *H 
pulses were optimized for 20 kHz of total offset range with RF miscalibration of ±5%, but the resulting pulses are robust for 
±15% RF miscalibration. *H PP 90° and UR 90° show the excitation profile for the z -> -y and 1 H PP 180° and UR 180° 
show an inversion profile from z — > — z. 
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Figure 16: The amplitude and phase of 13 C pulses are plotted as a function of the pulse duration. All are 1 ms long with a 
maximum amplitude of 10 kHz. 
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Figure 17: The experimental performance of 13 C pulses is plotted as function of the offset iw at ideal RF amplitude with a 
SF max of 10 kHz. Experiments were aquired on residual proton in HDO for offset range of 37 kHz in steps of 200 Hz. J H 
pulses were optimized for 35 kHz of total offset range with RF miscalibration of ±5% but the resulting pulses are robust for 
±10% RF miscalibration. 13 C PP 90° and UR 90° show the excitation profile for the z -> -y and 13 C PP 180° and UR 180° 
show an inversion profile from z — > — z. 
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Figure 18: The amplitude and phase of ^ and 13 C pulses optimized simultaneously for the magnetization transfers 2I Z S X — > 
2I y S x and 2I z S y —> — 2I y S z are plotted as a function of the pulse duration. Both pulses are 1 ms long with a maximum 
amplitude of 10 kHz. 
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13 offset points in each direction 

Figure 19: Simulated transfer efficiencies (TE) of the 1 H and 13 C pulses of Fig, |l8| optimized simultaneously for the magneti- 
zation transfers 2I Z S X — > 2I y S x (top) and 2I z S y — > —2I y S z (bottom), are plotted as a function of iff and 13 C offsets. 



28 



